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FOREWORD 


This  work  at  Massachusetts  Institute  of 
Technology  was  sponsored  by  the  Naval  Air 
Systems  Command  under  jnt.ract  No. 

N00019- 71-C-0220 .  The  contract  was  monitored 
by  Mr.  Ray  Malatino. 


ABSTRACT 


Computer  programs  have  been  developed  for  the  calcula¬ 
tion  of  helicopter  rotor  tip  vortx  geometry  in  hover  and 
forward  flight  and  for  the  calculation  of  helicopter  rotor 
harmonic  airloads  in  forward  flight.  Calculated  forward 
flight  tip  vortex  geometries  compare  well  in  general  with 
experimental  smoke  studies  although  there  are  differences 
in  detail.  The  hovering  tip  vortex  geometry  agrees  quali¬ 
tatively  with  experiment  but  does  not  move  downward  fast 
enough.  Airloads  were  computed  using  both  the  classical 
rigid  wake  assumption  and  the  distorted  tip  vortex  geometry 
obtained  from  the  computer  program.  When  compared  with 
experimental  airloads  measurements  the  rigid  wake  airloads 
give  better  results  than  the  distorted  wake  airloads.  This  is 
due  to  excessive  peak  heights  in  the  distorted  wake  airloads. 
The  problem  is  thought  to  be  the  neglect  of  some  unknown  but 
important  effects  which  limit  peak  heights  due  to  close  blade 
vortex  interactions. 

A  tv/o  dimensional  laser  Doppler  velocimeter  (LDV)  has 
been  developed  for  use  in  helicopter  rotor  wakes.  Extensive 
measurements  were  made  in  vortex  rings.  From  this  data  the 
circulation,  streamlines  and  vorticity  distribution  were  cal¬ 
culated.  Finally,  the  construction  of  a  hollow  composite 
material  rotor  for  use  with  the  laser  system  is  discussed. 
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SECTION  I 


HELICOPTER  ROTOR  WAKE  GEOMETRY  AND  AIRLOADS  WORK 

This  section  summarizes  research  done  on  the  computation 
jf  helicopter  rotor  wake  geometry  in  both  hovering  and  forward 
flight  and  the  effects  of  wake  geometry  on  rotor  harmonic 
airloads  in  forward  flight.  This  work  is  part  of  a  continuing 
program  of  research  in  this  area,  and  a  more  detailed  dis¬ 
cussion  of  the  results  of  this  work  plus  research  done  in 
previous  years  will  be  published  in  Reference  1. 

FORWARD  FLIGHT  MAKE  GEOMETRY 

The  only  concentrated  vortex  line  in  the  rotor  wake  is 
the  tip  vortex.  The  primary  source  of  high  harmonic  airloads 
and  hence  vibration  is  close  interactions  between  this  con¬ 
centrated  tip  vortex  and  the  rotor  blades.  In  order  to  pre¬ 
dict  the  airloads  due  to  these  close  interactions  the  geometry 
of  the  tip  vortex  must  be  known  accurately.  The  classical 
rigid  wake  geometry  assumption  based  on  uniform  inflow  is  not 
accurate  enough  for  this  purpose.  It  is  necessary  to  include 
the  distortion  of  the  tip  vortex  geometry  by  the  actual  non- 
uniform  and  time-varying  induced  velocities  experienced  by  the 
tip  vortex.  These  induced  velocities  are  a  function,  primarily, 
of  the  tip  vortex  geometry  which  is  in  turn  a  function  of  the 
induced  velocities.  Thus  the  computation  of  the  tip  vortex 
distortion  and  hence  the  geometry  require  the  integration  of 
the  Biot-Savart  relation  for  induced  velocities  over  the  wake 
which  must  be  repeated  (or  updated)  as  the  knowledge  of  the 
actual  geometry  of  the  wake  improves.  This  is  a  very  time- 
consuming  process  requiring  the  use  of  a  large  and  expensive 
computer  program. 

Digital  computer  programs  have  been  developed  at  M.I.T. 
and  elsewhere  to  compute  the  tip  vortex  distortion.  These 
programs  are  generally  more  than  an  order  of  magnitude  more 
expensive  than  the  corresponding  airloads  computation  for  a 
given  wake  geometry.  This  is  prohibitively  expensive  for  use 
in  routine  airloads  computations.  Research  effort  has  there¬ 
fore  been  concentrated  on  reducing  the  expense  of  wake  geometry 
computations  as  well  as  improving  their  accuracy. 

Most  of  the  expense  of  a  wake  geometry  computation  comes 
from  the  repeated  integrations  of  the  Biot-Savart  relation  for 
induced  velocities  over  the  wake.  Thus  the  best  way  to  save 
computation  is  to  maximize  the  efficiency  of  the  Biot-Savart 
relation  integrations  and  to  minimize  the  number  of  them.  The 
most  efficient  way  to  integrate  the  Biot-Savart  relation  over 
a  distorted  vortex  line  is  to  break  the  line  up  into  a  series 
of  short  straight  line  segments.  The  integration  can  be  performed 
in  closed  form  for  a  straight  line  segment.  Thus  the  integral 
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for  the  entire  line  is  obtained  by  adding  up  the  contributions 
of  all  of  the  straight  line  segments.  The  length  of  the 
straight  line  segments  is  a  trade-off  between  accuracy  and 
expense.  A  highly  distorted  tip  vortex  may  require  segments 
subtending  15°  or  less  in  azimuth  for  adequate  representation. 

To  minimize  the  number  of  straight  vortex  line  segments 
whose  contribution  to  the  induced  velocity  must  be  computed 
the  simplest  possible  wake  model  is  used.  The  wake  is  divided 
up  into  the  trailing  wake  due  to  radial  circulation  variations 
eind  the  shed  wake  due  to  azimuthal  circulation  variations. 

The  outboard  portion  of  the  trailing  wake  rolls  up  to  form  a 
concentrated  tip  vortex  which  is  the  dominant  feature  of  the 
wake  and  must  be  accurately  represented.  The  tip  vortex  is 
modeled  by  a  series  of  straight  vortex  line  segments  forming 
a  single  vortex  line.  These  vortex  line  segments  have  a 
finite  radius  viscous  core.  The  initial  core  radius  is  esti¬ 
mated  at  about  10%  of  the  blade  chord  based  on  experimental 
data.  After  a  segment  of  the  tip  vortex  undergoes  a  close 
interaction  with  a  rotor  blade  it  is  assumed  to  have  burst, 
increasing  the  core  radius  to  about  four  times  the  original 
radius . 

The  shed  wake  and  the  inboard  portion  of  the  trailing 
wake  remain  spread  out  in  vortex  sheets  near  the  rotor.  Due 
to  their  spread-out  nature  the  exact  location  cf  these  por¬ 
tions  of  the  wake  is  not  important  and  the  rigid  wake  geo¬ 
metry  assumption  can  be  used.  It  is  important,  however,  to 
avoid  introducing  unrealistic  induced  velocity  peaks  due  to 
the  use  of  powerful  concentrated  vortex  lines  to  represent 
these  portions  of  the  wake.  A  common  method  is  to  use  a 
large  number  of  vortex  lines  to  represent  the  inboard  trailing 
wake.  This  is  very  expensive  so  a  better  way  is  needed. 

For  the  most  accurate  representation  of  the  inboard  trail¬ 
ing  and  shed  wakes,  vortex  sheets  should  be  used.  The  Biot- 
Savart  relation  can  be  integrated  over  a  rectangular,  planar 
vortex  sheet  segment  in  closed  form.  The  inboard  trailing  and 
shed  wakes  can  therefore  be  represented  by  a  series  of  such 
vortex  sheet  segments.  The  first  attempts  to  use  vortex  sheets 
involved  representing  the  inboard  trailing  wake  with  one  set 
of  vortex  sheet  segments  and  the  shed  wake  with  two  sets,  one 
having  constant  radial  circulation  and  the  other  linear  radial 
circulation.  This  adds  up  to  three  vortex  sheet  segments  and 
one  vortex  line  segment  (for  the  tip  vortex) ,  per  azimuth  com¬ 
pared  to  three  vortex  line  segments  per  azimuth  for  the  all 
vortex  line  model.  Since  the  computation  of  the  induced  velo¬ 
city  contribution  of  a  vortex  sheet  segment  is  more  than  twice 
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as  expensive  as  the  same  computation  for  a  vortex  line  seg¬ 
ment  this  model  proved  to  be  nearly  three  times  as  expensive 
as  the  vortex  line  model.  The  results  of  tip  vortex  geometry 
computations  using  this  vortex  sheet  model  gave  better  agree¬ 
ment  with  experimental  results  than  computations  using  the 
vortex  line  model  however.  Since  the  circulation  of  the  shed 
wake  is  normally  almost  an  order  of  magnitude  less  than  the 
trailing  wake  a  model  was  tried  using  a  single  large  core 
vortex  line  with  linear  radial  circulation  for  the  shed  wake, 
plus  one  set  of  vortex  sheets  for  the  inboard  trailing  wake 
and  a  vortex  line  for  the  tip  vortex.  This  gives  the  same 
results  as  the  original  vortex  sheet  model  and  costs  only  40 
percent  more  than  the  vortex  line  model  and  has  been  adopted 
as  standard. 

After  choosing  a  wake  model  the  number  of  updatings  of  the 
integration  of  the  Biot-Savart  relation  must  be  minimized  to 
save  computation.  A  big  saving  can  be  made  by  dividing  the 
wake  elements  contributing  to  the  induced  velocity  at  a  point 
P  into  near  and  far  wake  elements  relative  to  point  P.  The 
induced  velocity  contribution  of  each  wake  element  varies  roughly 
inversely  with  the  square  of  the  distance  from  each  element  to 
point  P.  Therefore  the  induced  velocity  contributions  of  far 
wake  elements  at  point  P  do  not  need  to  be  updated  neaily  as 
often  as  the  contributions  of  the  near  wake  elements. 

The  M.I.T.  wake  geometry  program  is  iterative  starting  with 
the  rigid  wake  assumption.  An  iteration  ends  when  the  distor¬ 
tion  of  2vm  (typically  m  =  2)  in  azimuth  of  newly  generated 
wake  has  been  computed.  The  next  iteration  starts  assuming  the 
distortion  computed  by  the  previous  iteration.  The  induced 
velocity  contributions  of  each  far  wake  element  at  each  point  P 
are  computed  only  once  per  iteration.  The  induced  velocity 
contributions  of  the  near  wake  elements  are  updated  several 
times  during  an  iteration  but  not  at  every  azimuth  (time) 
interval. 

An  iteration  in  the  wake  geometry  computation  proceeds  by 
advancing  the  rotor  in  discrete  azimuth  steps  thus  generating 
new  wake.  The  induced  velocity  and  hence  the  tip  vortex  dis¬ 
tortion  is  computed  only  at  the  boundary  between  the  new  wake 
generated  during  the  current  iteration  and  the  old  wake  from 
the  previous  iteration.  As  more  wake  is  generated  its  contri¬ 
bution  to  the  induced  velocity  is  added  in  while  the  contribu¬ 
tions  of  the  rest  of  the  wake  are  assumed  to  remain  constant. 
After  n  azimuth  steps  the  induced  velocity  contributions 
from  the  entire  wake  are  recomputed  (updated)  to  account  for 
the  movement  of  the  boundary  point  under  the  influence  of  the 
previously  computed  induced  velocities  and  to  account  for  the 
improved  knowledge  of  the  wake  geometry.  This  updating  is  done 
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independently  for  the  boundary  point  at  each  different  azimuth. 
The  updating  parameter  n  is  a  function  of  azimuth  and  smaller 
values  are  used  at  the  sides  of  the  rotor  where  the  tip  vortices 
of  the  various  blades  and  successive  turns  of  each  tip  vortex 
cross/  resulting  in  large  tip  vortex  distortions.  In  addition 
the  tip  vortex  distortion  of  all  of  the  wake  generated  so  far 
in  the  current  iteration  is  updated  every  l  azimuth  steps. 

This  involves  recomputing  the  induced  velocity  and  hence  the 
distortion  at  all  of  the  points  in  the  new  wake,  not  just  at 
the  boundary  between  the  old  and  the  new  wake.  This  is  much 
more  expensive  than  simply  updating  the  induced  velocity  at 
the  boundary  and  hence  is  not  done  as  frequently.  Since  the 
entire  new  wake  is  involved,  this  type  of  updating  is  done  at 
all  azimuths,  so  Z  is  not  a  function  of  azimuth. 

To  speed  up  convergence  of  the  wake  geometry  it  is  impor¬ 
tant  to  use  the  best  available  knowledge  of  the  wake  geometry 
at  all  points  in  the  computation.  Thus,  if  the  contribution 
of  a  wake  element  in  the  old  wake  is  being  computed,  the  dis¬ 
tortion  at  that  wake  element  is  determined  by  taking  the 
distortion  accumulated  at  that  azimuth  up  to  its  arrival  at 
the  boundary  between  the  old  and  new  wakes, as  computed  for  the 
new  wake, and  adding  in  the  old  wake  distortion  for  the  time  after 
it  left  the  boundary.  This  same  principle  is  applied  in  the 
middle  of  an  update  when  the  first  part  of  the  new  wake  has 
been  updated  but  the  last  part  has  not.  Since  the  old  wake 
does  not  extend  to  infinity,  and  in  fact  is  just  as  long  as 
the  new  wake  at  the  end  of  an  iteration,  extrapolation  is  used 
to  find  the  geometry  of  wake  elements  beyond  the  end  of  the 
old  wake.  This  ensures  that  there  is  always  a  realistic  repre¬ 
sentation  of  the  wake  beyond  the  boundary  between  the  old  and 
new  wakes,  even  when  this  boundary  reaches  the  lower  end  of  the 
old  wake. 

During  wake  geometry  computations  using  this  program  the 
location  of  some  wake  elements  tends  to  oscillate.  These 
oscillations  may  take  six  or  more  iterations  to  converge  to  a 
final  solution,  which  is  very  expensive.  To  damp  these  oscil¬ 
lations  more  rapidly  the  first  step  is  to  average  the  distor¬ 
tions  obtained  from  successive  iterations.  This  is  still  too 
slow,  so  the  averaging  of  the  distortions  obtained  from 
successive  updates  is  used.  This  normally  results  in  conver¬ 
gence  within  two  iterations. 

The  wake  is  not  the  only  source  of  induced  velocity  and 
hence  distortion.  The  bound  vorticity  also  makes  a  signifi¬ 
cant  contribution.  The  bound  vorticity  is  represented  by  a 
single  vortex  line  at  the  quarter  chord  line  of  the  blade. 

The  distortion  is  computed  at  points  every  15  degrees  or  so 
along  the  tip  vortex.  In  cases  where  the  tip  vortex  passes 
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close  under  or  over  the  blade  and  roughly  perpendicular  to  it 
the  contribution  of  the  bound  vorticity  to  the  distortion  is 
largely  determined  by  the  distance  from  the  blade-tip  vortex 
crossing  to  the  nearest  point  on  the  tip  vortex  at  which  the 
distortion  is  computed.  Since  these  points  are  about  15  degrees 
or  a  quarter  of  the  rotor  radius  apart,  this  results  in  large, 
random  variations  in  the  contribution  of  the  bound  vorticity  to 
the  distortion.  The  way  to  avoid  this  problem  is  to  actually 
integrate  over  time  the  induced  velocity  at  a  point  P  on  the 
tip  vortex  as  the  blade  passes  under  or  over  it,  instead  of 
computing  the  induced  velocity  at  two  different  times  15  degrees 
in  azimuth  apart  and  multiplying  the  average  by  delta  time.  It 
turns  out  that  representing  the  bound  vorticity  by  a  vortex 
sheet  segment  15  degrees  wide  is  the  same  as  performing  this 
integral  directly.  Therefore  the  vortex  sheet  representation 
of  the  bound  vorticity  is  used  whenever  the  tip  vortex  passes 
close  over  or  under  a  blade.  This  should  not  be  confused  with 
a  lifting  surface  representation  of  the  bound  vorticity.  It  is 
only  a  method  of  integrating  the  effect  of  a  lifting  line  over 
t  ime . 


The  wake  geometry  program  has  been  tested  for  1,  2  and  4 
bladed  cases.  The  results  shov.  better  agreement  with  experi¬ 
ment  than  the  old  wake  geometry  program  developed  under  the 
previous  contract.  In  addition  faster  convergence  and  faster 
execution  have  been  achieved.  The  developments  which  produced 
these  improvements  are:  the  modified  vortex  sheet  wake  model, 
the  inclusion  of  vortex  core  bursting,  the  vortex  sheet  repre¬ 
sentation  of  the  bound  vorticity,  improved  updating  procedures 
including  making  n  a  function  of  azimuth,  improved  procedures 
for  always  using  the  most  recent  wake  geometry  information, 
and  averaging  the  distortion  from  successive  updates. 

The  wake  geometry  program  has  been  run  for  various  values 
of  the  computational  parameters  n,  £,  and  d^  (a  parameter  which 

determines  the  split  between  the  near  and  far  wakes).  Values 
of  these  parameters  have  been  found  which  give  a  good  cost- 
accuracy  trade-off.  In  addition  the  effects  of  vortex  core 
size  and  various  wake  models  have  been  evaluated.  These 
results  will  be  presented  in  more  detail  in  Reference  1. 

EFFECTS  OF  WAKE  GEOMETRY  ON  HARMONIC  AIRLOADS 


The  primary  reason  for  computing  the  tip  vortex  distortion 
is  to  improve  the  accuracy  of  harmonic  airloads  computations. 
Therefore  airloads  computations  have  been  made  using  both  rigid 
and  distorted  wake  models  for  2,  4  and  6  bladed  rotors  for  which 
experimental  airloads  data  are  available.  Unfortunately  there 
is  no  case  for  which  both  experimental  airloads  and  tip  vortex 
geometry  are  available.  However,  comparison  with  experimental 
tip  vortex  geometry  data  for  cases  where  no  experimental  airloads 
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data  are  available  has  shown  that  the  results  of  the  wake 
geometry  program  are  generally  in  agreement  although  there  are 
local  differences. 

The  wake  geometry  program  is  normally  run  in  conjunction 
with  the  airloads  program.  The  first  step  is  to  make  a  rigid 
wake  airloads  run.  The  circulation  distribution  and  the  loca¬ 
tion  of  any  vortex  core  bursting  is  input  from  this  run  into 
the  wake  geometry  program.  The  tip  vortex  geometry  obtained 
from  this  run  is  then  used  in  a  second  airloads  run.  This 
process  continues  as  long  as  necessary.  Normally  one  wake 
geometry  run  and  two  airloads  runs  are  sufficient,  but  some¬ 
times  the  circulation  distribution  and/or  the  distribution 
of  vortex  core  bursting  may  change  radically  due  to  tip 
vortex  distortion.  In  this  case  a  second  wake  geometry  run 
and  a  third  airloads  run  will  be  needed. 

The  airloads  program  uses  a  wake  model  similar  to  the 
large  core  vortex  line  model  of  the  wake  geometry  program 
including  vortex  core  bursting.  Under  the  next  contract 
the  airloads  program  will  be  modified  to  use  a  single  vortex 
sheet  plus  two  vortex  lines  wake  model  similar  to  the  latest 
wake  geometry  program.  A  vortex  sheet  segment  is  used  for  the 
shed  wake  element  immediately  behind  the  blade  and  the  inte¬ 
gration  over  the  shed  wake  is  carried  to  within  e  of  the 
control  points  on  the  blade,  where  e  is  determined  from 
lifting  surface  theory. 

Whenever  the  tip  vortex  passes  close  under  or  over  a 
blade  a  rapid  radial  variation  in  lift  and  hence  circulation 
is  produced  on  the  blade.  This  in  turn  results  an  extra  trail¬ 
ing  wake  being  generated.  This  extra  trailing  wake  could  be 
represented  by  a  many  ( ~20 ) trailing  vortex  line  near  wake  model 
but  this  is  expensive  and  neglects  lifting  surface  effects 
which  become  important  for  very  close  blade-tip  vortex  inter¬ 
actions.  Instead  the  approximations  developed  in  Reference  2 
using  lifting  surface  theory  are  used  whenever  the  tip  vortex 
passes  close  over  or  under  a  blade.  In  addition  the  vortex 
core  is  assumed  to  burst  (increase  by  about  a  factor  of  4  in 
size)  after  a  close  blade-tip  vortex  interaction. 

Once  the  induced  velocity  field  has  been  determined  the 
airloads  and  circulation  distribution  are  computed.  Rigid  body 
flapping  and  first  mode  flapwise  bending  are  included  but  lag 
and  torsion  are  neglected.  Provision  is  made  for  flapping 
hinge  offset  and  a  spring  at  the  flapping  hinge.  Reverse 
flow,  radial  flow,  and  dynamic  stall  effects  are  included. 
Iteration  is  used  to  obtain  a  converged  circulation  distribu¬ 
tion.  To  save  computation  a  matrix  of  influence  coefficients 
is  computed  for  a  given  wake  geometry,  which,  when  multiplied 
by  the  circulation  matrix,  gives  the  induced  velocity  matrix. 
Most  of  the  computation  in  an  airloads  run  is  in  finding  the 
influence  coefficient  matrix.  The  circulation  iteration 
using  that  matrix  is  very  fast. 
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The  case  most  carefully  studied  so  far  is  a  four  bladed 
rotor  at  an  advance  ratio  of  0.18.  For  this  case  the  rigid 
wake  airloads  computation  gives  good  overall  results  when 
compared  with  experimental  data.  The  distorted  wake  airloads 
computation  improves  agreement  in  some  local  areas  but  overall 
is  much  worse  than  the  rigid  wake.  The  basic  problem  is 
larger  peaks  in  the  theoretical  airloads  than  in  the  experimental 
ones.  The  distorted  wake  is  worse  than  the  rigid  wake  because 
the  wake  is  closer  to  the  rotor,  especially  over  the  forward 
edge  of  the  rotor.  The  tip  vortex  geometry  is  probably  generally 
correct,  with  possible  local  exceptions,  because  the  tip  vortex 
should  be  closer  to  the  rotor  than  rigid  wake  would  predict 
over  the  forward  edge  of  the  rotor,  since  this  is  an  area  of 
upwash  from  the  tip  vortices  of  the  preceding  blades.  There¬ 
fore  the  problem  is  inadequate  treatment  of  close  blade-tip 
vortex  interactions,  despite  all  of  the  special  provisions 
made  for  this  so  far. 

The  better  results  obtained  by  the  rigid  wake  model  are 
due  to  compensating  errors,  the  reduction  in  peak  size  caused 
by  the  wake  being  too  far  away  from  the  blades  being  compen¬ 
sated  by  neglect  of  close  blade-tip  vortex  interaction  effects 
which  if  included  would  have  reduced  peak  sizes.  The  most 
likely  close  blade-tip  vortex  effects  are:  local  separation 
or  stall  on  the  blade  due  to  the  intense  induced  velocities 
resulting  in  steep  local  pressure  gradients,  vortex  burst¬ 
ing  over  or  under  the  blade  instead  of  behind  it,  the  possibility 
that  seme  of  the  tip  vortex  circulation  is  not  wrapped  up  into 
a  tight  vortex  core  but  is  spread  out  in  spiral  vortex  sheets 
around  the  core,  and  local  tip  vortex  distortion  near  the  blade. 
The  first  two  effects  are  being  studied  in  a  flow  visualiza¬ 
tion  program  at  M.I.T.  The  third  effect  is  suggested  by  the 
fact  that  some  investigators  find  only  60%  or  so  of  the  bound 
circulation  in  the  tip  vortex  core,  for  example,  Reference  3. 

The  wake  geometry  program  is  only  intended  to  provide  the 
gross  tip  vortex  geometry.  It  only  computes  the  tip  vortex 
distortion  at  points  15  degrees  or  so  apart  and  a  lifting  line 
is  used  to  represent  the  blade.  To  compute  the  local  tip  vor¬ 
tex  distortion  near  a  blade  the  rest  of  the  wake  can  be  neglected 
but  the  blade  must  be  more  carefully  represented.  At  least  a 
lifting  surface  representation  will  be  needed  and  thickness 
effects  may  turn  out  to  be  important  as  well. 

In  summary,  the  rigid  wake  model  gives  good  results  due  to 
compensating  errors.  To  improve  on  the  rigid  wake  results 
it  will  be  necessary  to  use  distorted  tip  vortex  geometry  in 
conjunction  with  a  much  more  careful  treatment  of  close  blade- 
tip  vortex  interaction  effects. 


HOVERING  WAKE  GEOMETRY 


A  specialized  version  of  the  wake  geometry  program  has 
been  developed  for  the  hovering  case.  This  program  uses  the 
assumption  of  azimuthal  symmetry  to  save  computation.  Since 
there  is  no  inflow  due  to  forward  speed,  and  since  the  wake  is 
not  swept  downstream  by  the  forward  speed,  more  of  the  wake 
is  important  than  for  a  forward  flight  wake  geometry  compu¬ 
tation.  Therefore  an  infinite  vortex  cylinder  is  added  to 
the  wake  model  starting  at  the  lower  end  of  the  normal  wake 
model  and  continuing  to  infinity.  Due  to  the  assumption  of 
azimuthal  symmetry  there  is  no  shed  wake.  In  addition  to  the 
standard  options  of  a  rigid  wake  inboard  trailing  vortex  line 
or  vortex  sheet,  as  in  forward  flight,  there  is  provision  for 
a  tilted  inboard  trailing  vortex  sheet  with  its  outer  edge 
moving  down  faster  than  its  inner  edge.  This  last  option  is 
based  on  experimental  data  on  the  inboard  trailing  sheet  geo¬ 
metry  from  Reference  4.  Other  than  the  above  mentioned  modifi¬ 
cations  the  hovering  wake  geometry  program  is  similar  to  the 
forward  flight  wake  geometry  program  described  previously. 

Great  difficulty  was  experienced  in  obtaining  converged 
solutions  from  the  hovering  wake  geometry  program.  Averaging 
of  the  distortion  from  successive  updates,  as  described  pre¬ 
viously  for  the  forward  flight  program,  was  originally 
developed  for  the  hovering  case.  Using  this  averaging,  plus 
very  careful  treatment  of  the  wake  beyond  the  lower  bound  of 
wake  geometry  computation,  including  the  infinite  vortex 
cylinder,  fairly  good  convergence  is  obtained  after  6  to  8 
iterations.  The  resulting  wake  geometry  does  not  move  down  as 
fast  or  contract  as  fast  as  the  experimental  results  of 
Reference  4  indicate  it  should. 

In  addition  some  recent  unpublished  theory  by  Professor 
Widnall  at  M.I.T.  indicates  that  an  infinite  helical  vortex 
cylinder  is  always  unstable.  Since,  after  the  first  turn  or 
so,  the  hovering  tip  vortex  is  a  semi-infinite  helical  vortex 
cylinder,  it  is  also  likely  to  be  unstable.  This  explains 
the  difficulty  in  obtaining  a  converged  solution.  The  con¬ 
verged  solutions  that  were  finally  obtained  are  based  on  the 
use  of  averaging  of  the  distortion  between  successive  updates. 
This  averaging  is  in  turn  based  on  the  assumption  that  a  con¬ 
verged  solution  exists.  The  net  result  is  to  introduce  some 
extra  damping  into  the  system  which  does  not  exist  in  the  real 
world.  In  addition  the  wake  geometry  is  only  computed  every 
20  degrees  in  azimuth  which  means  that  short  wave  length 
instabilities  do  not  show  up.  The  converged  tip  vortex  geo¬ 
metry  produced  by  this  program  probably  represents  some  sort 
of  average  location  of  the  tip  vortex  and  may  be  useful  as  such. 


The  tip  vortex  near  the  rotor  is  the  part  of  the  wake  which 
is  important  for  airloads  or  performance  computations.  This 
part  of  the  wake  appears  to  be  stable  in  the  real  world  as 
shown  in  the  experimental  results  of  Reference  4.  Therefore 
effort  should  be  concentrated  on  computing  the  geometry  of  this 
part  of  the  tip  vortex.  The  results  of  the  current  hovering 
wake  geometry  program  do  not  go  down  fast  enough  or  contract 
fast  enough  in  this  area,  when  compared  to  experiment.  Hand 
computations  show  that  this  problem  is  very  sensitive  to  small 
changes  in  the  geometry,  since  the  net  induced  velocity  on  an 
element  of  the  tip  vortex  is  normally  the  small  difference  of 
large  numbers.  This  suggests  that  even  a  very  small  amount  of 
wind  or  recirculation  could  have  significant  effects  on  the  tip 
vortex  geometry. 

In  summary,  the  current  hovering  wake  geometry  program 
developed  from  the  forward  flight  wake  geometry  program  does 
not  give  satisfactory  results.  A  completely  new  hovering  wake 
geometry  program  should  be  developed,  independent  of  the  for¬ 
ward  flight  program,  which  concentrates  on  computing  the 
geometry  of  the  first  part  of  the  tip  vortex.  This  new  pro¬ 
gram  should  be  designed  to  separate  our.  the  effects  of  the 
various  contributors  to  tip  vortex  distortion  (i.e.  the  bound 
vorticity,  the  curvature  of  the  tip  vortex,  the  tip  vortices 
of  the  other  blades,  the  returning  part  of  the  tip  vortex 
itself,  the  inboard  trailing  wake,  etc.).  Using  this  tool,  a 
better  understanding  of  the  relative  importance  of  the  various 
parts  of  the  problem  can  be  obtained  and  used  to  develop  a 
solution. 


REFERENCES 


Scully,  M.  P. ,  "Computation  of  Helicopter  Rotor  Wake  Geometry 
and  Its  Influence  on  Rotor  Harmonic  Airloads,"  M.I.T.  Dept 
of  Aeronautics  and  Astronautics,  PhD  Thesis  (to  be  published). 

Johnson,  W.  R. ,  "A  Lifting  Surface  Solution  for  Vortex 
Induced  Airloads  and  Its  Application  to  Rotary  Wing  Airloads 
Calculations,"  M.I.T.  Aeroelastic  and  Structures  Research 
Laboratory,  Technical  Report  TR  153-2,  April  1970. 

Chigier,  N.  A.  and  Corsiglia,  V.  R. ,  "Tip  Vortices  - 
Velocity  Distributions,"  27th  Annual  National  V/STOL  Forum 
of  the  American  Helicopter  Society,  Preprint  S22,  May  1971. 

Landgrebe,  Anton  T.  ,  "An  Analytical  and  Experimental 
Investigation  of  Helicopter  Rotor  Hover  Performance  and 
Wake  Geometry  Characteristics,"  USAAAMRDL  Tech.  Report 
71-24,  June  1971. 


10 


SECTION  II 


DEVELOPMENT  OF  LASER  DOPPLER  VELOCIMETER 
FOR  USE  IN  HELICOPTER  ROTOR  WAKES 


INTRODUCTION 


Traditionally  the  velocity  of  a  time  varying  flow  field 
has  been  measured  using  a  hot-wire  anemometer.  The  use  of  a 
hot  wire  in  the  wake  of  a  helicopter  (or  any  flow  containing 
concentrated  vorticity)  poses  two  serious  problems:  probe 
interference  (possibly  in  the  form  of  vortex  breakdown)  and 
calibration  in  a  flow  field  whose  direction  and  magnitude  is 
rapidly  varying.  For  these  reasons  a  laser  Doppler  veloci- 
meter  (LDV)  has  been  developed. 


First  the  basic  principles  involved  in  an  LDV  system  are 
discussed  and  then  a  two-dimensional  LDV  system  is  described. 
This  system  is  used  for  extensive  measurements  on  vortex  rings 
Finally  the  construction  of  a  hollow  composite  material  rotor 
for  use  with  the  LDV  system  is  discussed. 


BASIC  PRINCIPLES  OF  LASER  DOPPLER  VELOCIMETRY 


Since  the  introduction  of  the  laser  Doppler  velocimeter 
(LDV)  in  1965  by  Yeh  and  Cummings  (Ref  1)  numerous  systems 
have  been  reported.  Two  of  thesQ,  the  cross  beam  —  reference 
beam  system  (Fig.  la)  of  ref  2,  and  the  dual  scatter  system 
(Fig.  lb)  of  refs  3  and  4  will  be  analyzed.  These  two  were 
chosen  because  they  are  self  aligning;  that  is,  once  the  beams 
are  made  to  intersect  at  the  probing  point  and  pass  throug. 
the  aperature  in  front  of  the  photomultiplier  tube,  the  system 
is  aligned.  This  eliminates  the  tedious  alignment  and  vibration 
problems  of  other  systems. 


Laser  light  that  is  scattered  from  particles  (smoke,  dirt, 
polystyrene  spheres)  moving  at  th  local  fluid  velocity*  under¬ 
goes  a  Doppler  shift  given  by 


f  =  4  o(f!) 


The  difference  between  the  fluid  and  particle  velocities  is 
estimated  in  the  Appendix. 


r\ 


where 


f  /  X,  kj  =  frequency,  wavelength  and  unit  vector  in 
1  the  propagation  direction  of  the  incident 
laser  beam 

V  =  velocity  of  particle 

A 

kg  =  unit  vector  in  the  direction  that  the 

scattered  light  is  detected 

c  =  velocity  of  light 

The  frequency  of  the  scattered  radiation  f,  which  is  of  the 
order  of  light  frequencies  (^5  x  10^  cy/sec)  can  be  detected 
either  directly  using  a  Fabry-Perot  interferometer  or  can  be 
heterodyned  with  some  light  wave  at  a  known  frequency  to 
obtain  a  signal  whose  frequency  is  equal  to  the  shift  of  the 
laser  light  frequency.  The  present  resolution  of  Fabry-Perot 
interferometers  is  only  ^2MHz  which  implies  a  velocity  reso¬ 
lution  of  ^10  ft/sec.  For  the  present  study  this  is  inadequate 
so  a  heterodyning  system  is  used. 

In  the  two  types  of  heterodyne  systems  studied  the  cross¬ 
beam  reference  beam  system  (Fig.  la)  heterodynes  the  scattered 
radiation  with  radiation  at  the  xaser  frequency  fQ  and  the 

dual  scatter  system  (Fig.  lb)  heterodynes  the  scattered  radia¬ 
tion  of  two  different  incident  beams.  In  both  cases  the 
difference  frequency  is  calculated, using  Eq.  1  and  the  geometry 
of  Fig.  1,  to  be 


c  _  2  U  sin  0/2 
D  X 


where  U  is  the  velocity  component  in  a  direction  perpendicular 
to  the  bisector  of  the  angle  between  two  beams.  Fig. 2  0 
shows  fQ  versus  U  for  various  values  of  0  with  a  =  6328A,  the 
wavelength  of  a  He  Ne  laser.  It  should  be  noted  that  there  is 
a  180°  ambiguity  in  determining  direction  of  U.  The  plus  or 
minus  direction  of  U  must  be  determined  by  prior  knowledge  or 
flow  visualization. 


The  output  signal  of  the 
proportioned  to  the  intensity 
system  is 


photomultiplier  (P.M.)  tube  is 
which  for  the  reference  beam 


I 


Re  = 


v  Ep/  E*(t)  cos(z-tr^t') 
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when  a  single  particle  traverses  the  measuring  volume.  This 

produces  an  output  like  that  shown  in  Fig.  3a.  The  d.c.  level 

is  essentially  constant  since  E  (t)  < <  E  .  The  a.c.  component 

s  c 

is  modulated  by  the  term,  EE  (t)  which  is  generally  a  gaussian. 

K  S 

In  the  dual  scatter  system  the  light  scattered  from  one 
beam  is  heterodyned  with  light  scattered  from  the  other  beam 
resulting  in 


|E^r+  |E^f*  ES|(t)-ESi(t)co^TT+0t) 


The  results  for  E  (t)  =  E  (t)  are  shown  in  Fig.  3b.  The 

S1  s  2 

d.c.  component  is  now  replaced  by  a  low  frequency  component. 

The  actual  signal  when  a  .0002  wire  traverses  the  measuring 
volume  is  shown  in  the  oscilloscope  trace  of  Fig.  4a.  When 
the  wire  doesn't  pass  through  the  center  of  the  measuring 
volume  a  signal  (Fig.  4b)  containing  a  large  amount  of  the  low 
frequency  component  and  small  amount  of  the  Doppler  signal 
component  is  obtained. 

In  the  present  investigation  more  than  one  scatterer  will 
be  in  the  measuring  volume  at  one  time  resulting  in  a  signal 
which  is  the  sum  of  the  signals  from  individual  particles.  A 
typical  "continuous"  signal  for  a  dual  scatter  system  is  shown 
in  Fig.  4c.  This  signal  was  obtained  from  a  rotating  ground 
glass  disk.  It  shows  the  general  features  of  a  dual  scatter¬ 
ing;  namely,  a  low  frequency  component  and  low  frequency  modu¬ 
lation  of  the  Doppler  signal.  The  low  frequency  modulation  is 
a  result  of  random  constructive  and  destructive  interference 
among  the  waves  scattered  by  many  particles,  causing  the  sig¬ 
nal  to  go  randomly  to  zero  during  severe  destructive  interfer¬ 
ence.  This  phenomenon,  called  "drop-out",  and  the  low  fre¬ 
quency  component,  are  the  major  problems  in  the  electronic 
processing  of  the  signal. 

A  comparison  of  the  frequency  spectrums  of  the  two  systems 
is  shown  in  Fig.  5  for  scattering  from  a  ground  glass  disk  and 
a  turbulent  jet.  Although  the  magnitudes  of  the  signals  are 
about  equal,  the  low  frequency  components  in  the  dual  sc  itter 
system  is  greater. 

The  object  of  this  investigation  is  to  determine  which  of 
the  two  systems  will  give  a  better  signal-to-noise  ratio.  The 
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reference  beam  system  has  the  advantage  that  the  signal  can 
be  increased  by  increasing  the  strength  of  the  reference 
beam.  However,  the  solid  angle  through  which  scattered 


light  may  be  collected  is  limited  since  the  Doppler  shift 

— is  a  function  of  viewing  direction  k  so  that 
A-  s 

at  large  solid  angles  the  Doppler  frequency  would  be 

unacceptably  broadened.  On  the  other  hand,  the  frequency 

A 


of  a  dual  scatter  system  depends  only  on 


the  incident  beam  directions  kT,,kT  so  that  the  solid  angle 

1  2 


for  viewing  scattered  radiation  is  unlimited.  However,  the 
strength  of  the  signal  from  a  dual  scatter  system  is  deter¬ 
mined  entirely  by  the  scattering  t.ijcess.  Also  there  will  be 
more  low  frequency  components  in  the  dual  scatter  system. 


In  order  to  complete  the  calculations  the  quantities  Eg , 

E  and  E  must  be  found  which  requires  a  detailed  study  of 
S1  s2 

the  actual  scattering  process.  This  study  is  not  yet  complete 
but  will  be  reported  at  a  later  date. 


The  effect  of  polarization  on  the  intensity  of  the  signal 
is  seen  by  letting  the  direction  of  polarization  of  the 
incident  beams  be  given  by  the  unit  vector  e  with  appropriate 
subscript.  It  is  assumed  and  experiments  confirm  that  the 
polarization  is  not  changed  by  the  scattering  process,  so  for 
the  dual  scatter  system 

*  E  A 


where  the  directions  a,  b  are  shown  in  Section  AA  of  Fig.  1. 

The  Doppler  shift  term  of  the  intensity  is  proportional 
to 

ESi-Ek=  L"  E«,E.xcos©  +  LbiEbtj  (2) 

For  arbitrary  polarization,  Eq.  1  shows  that  the  signal 
decreases  as  0  increases.  If  the  incident  beams  are  linearly 

polarized  in  the  b.and  b„  (i.e.  E  =  E  =0)  directions 

1  2  ai  a2 

then  the  maximum  signal  will  be  obtained. 
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At  the  present  time,  the  investigation  does  not  indicate 
which  system  is  better.  The  experiments  indicate  that  both 
systems  perform  well. 

TWO  DIMENSIONAL  LDV  SYSTEM  AND  ASSOCIATED  ELECTRONICS 

The  LDV  system  constructed  for  extensive  measurements  is 
a  dual  scatter  system.  A  one-dimensional  dual  scatter  system 
is  shown  in  Fig.  6.  The  laser  beam  is  divided  by  a  standard 
prism-type  beamsplitter  with  a  non-absorbing  beamsplitting 
interface.  This  method  of  dividing  the  laser  beam  is  prefer¬ 
able  to  other  systems  using  mirrors  and  beamsplitters  since 
only  one  standard  optical  component  is  needed  and  the  path 
lengths  (distance  from  laser  to  measuring  point)  of  two  beams 
are  equal. 

A  two-dimensional  system  is  constructed  by  adding  a  second 
beamsplitter  (rotated  90  degrees  with  respect  to  the  first)  to 
form  four  parallel  beams.  The  pattern  of  the  four  beams  is 
shown  in  Fig.  7.  The  polarization  of  the  beams  is  chosen  to 
maximize  the  signal  (Eq.  1)  and  to  eliminate  crosstalk  between 
the  two  dimensions.  By  placing  an  analyzing  polaroid  in  front 
of  the  P.M.  tube  either  of  the  velocity  components  can  be  chosen 
with  no  interference  from  the  other  component.  This  is  possible 
because  the  polarization  is  changed  only  slightly  by  the  scatter¬ 
ing  process.  In  the  present  system,  both  components  are  pro¬ 
cessed  simultaneously  through  one  P.M.  tube  using  the  electronic 
system  described  below.  A  photograph  of  the  two-dimensional 
LDV  system  is  shown  in  Fig.  8.  The  two  components  could  be 
processed  separately  by  splitting  the  scattered  radiation  with 
a  polarizing  type  beamsplitter  and  using  two  P.M.  tubes  and  two 
sets  of  electronics. 

The  electronic  system  for  processing  the  LDV  signal  is  one 
suggested  by  Ref.  5  for  use  with  periodic  flow  fields.  This 
"sampling"  system  will  be  described  for  use  with  periodically 
produced  vortex  rings  in  terms  of  the  time  diagrams  of  Fig.  9 
and  the  equipment  layout  of  Fig.  6.  The  function  generator 
initiates  a  pulse  which  is  amplified  and  used  to  drive  a  loud¬ 
speaker  which  forms  a  vortex  ring  at  a  sharp  edged  orifice.  A 
synchronous  pulse  is  also  taken  off  of  the  function  generator 
(Fig.  9a) . 

The  velocity  versus  time,  as  the  vortex  ring  goes  past  the 
measuring  point,  is  shown  in  Fig.  9b.  This  is  proportional  to 
the  frequency  versus  time  of  the  output  signal  of  the-  P.M.  tube 
which  is  fed  into  the  input  of  the  spectrum  analyzer.  This 
synchronous  pulse  from  the  function  generator  initiates  a  time 
delay  which  starts  the  sawtooth  sweep  generator  (Fig.  9c)  for 
the  spectrum  analyzer  after  a  time  t^.  The  linear  sawtooth 
represents  the  center  frequency  of  uthe  narrow  filter  of  the 
spectrum  analyzer.  When  the  frequency  from  the  P.M.  tube  is 
equal  to  the  instantaneous  frequency  of  the  narrow  band  pass 
filter  of  the  spectrum  analyzer  a  "pip"  will  occur  on  the  spec¬ 
trum  analyzer  face  and  at  the  spectrum  analyzer  output  (Fig.  9d) 
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This  condition  can  be  determined  by  superposing  Fig.  9c  on  Fig. 

9b.  The  linear  sawtooth  is  fed  to  ,;he  input  of  a  gate  which 
is  only  opened  when  a  "pip"  occurs  on  the  spectrum  analyzer. 

The  output  of  the  gate  --  in  this  case,  three  spots  —  is  stored 
on  a  storage  oscilloscope.  The  three  spots  correspond  to  a  zero 
frequency  marker  and  two  points  on  the  velocity  versus  time 
curve  we  wish  to  determine.  By  varying  the  time  delay  the  re¬ 
maining  points  on  the  velocity  curve  can  be  stored. 

The  results  of  demodulating  a  test  oscillator  signal  vary¬ 
ing  from  50  kHz  to  950  kHz  at  a  rate  of  33  Hz  are  shown  in  Fig.  10, 

The  entire  LDV  system  was  checked  by  measuring  two  components 
of  the  velocity  of  a  spinning  disk  as  the  disk  was  traversed 
across  the  measuring  point.  The  x  component  of  velocity  U  =  -  wy 

should  be  a  constant  and  the  y  component  V  =  uix  should  vary 

linearly  as  the  disk  moves  in  the  x  direction.  The  results  of 

this  check,  shown  in  Fig.  11,  indicate  that  the  linearity  of  the 

electronics  system  is  excellent  and  that  the  velocity  error  is 
within  the  error  in  determining  the  r.p.m.  of  the  disk.  The 
ambiguity  in  determining  the  direction  of  U  and  V  is  demonstrated 
in  Fig.  H.  The  U  component  is  always  negative  and  the  V  compo¬ 
nent  is  negative  for  negative  values  of  X. 

The  two  dimensional  system  has  been  used  to  measure  the  velo¬ 
city  distribution  in  periodically  proaiced  vortex  rings  (smoke 
rings) .  (The  results  of  extensive  measurements  are  presented  in 
the  next  section.)  The  rings  are  similar  to  a  helicopter  rotor 
wake  since  regions  of  time  varying  concentrated  vorticity  are 
produced,  yet  the  flow  field  is  simpler  ^ince  it  is  axisymmetric . 

A  picture  of  a  vortex  ring  is  shown  in  Fig.  12  and  the  corres¬ 
ponding  velocity  distributions  in  Fig.  13.  The  velocity  curves 
are  the  U  and  V  components  of  velocity  versus  time  as  the  vortex 
ring  travels  by  the  fixed  measuring  point.  There  is  a  180° 
ambiguity  in  the  velocity  vector  (since  frequency  is  always  posi¬ 
tive)  .  Flow  visualization  indicates  that  the  U  component  is 
always  positive  while  the  V  component  is  positive  to  the  left  of 
the  central  zero  and  negative  to  the  right. 
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TWO  DIMENSIONAL  MEASUREMENTS  OF  VORTEX  RINGS 

Once  the  principle  of  using  an  LDV  system  for  measuring  velo¬ 
cities  in  vortex  rings  was  demonstrated,  detailed  surveys  of  two 
vortex  rings  (denoted  Ring  No.  1  and  Ring  No.  2)  were  made.  Com¬ 
plete  two  dimensional  data  were  taken  for  Ring  No.  1  but  only  the 
component  of  velocity  in  the  direction  of  ring  travel  was  measured 
for  Ring  No.  2. 

The  data  are  treated  in  the  following  quasi-steady  manner. 

The  LDV  system  measures  velocity  versus  time  at  a  fixed  point 
in  space  as  the  vortex  ring  goes  by  the  measuring  point.  The 
velocity  distribution  is  transformed  to  the  steady  coordinate 
system  of  Fig.  14  by  letting  Z  =  UQt  where  UQ  is  the  instantan¬ 
eous  velocity  of  the  vortex  ring  at  the  measuring  point.  The 
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quasi-steady  approximation  assumes  that  the  properties  of  the 
vortex  ring  do  not  change  (due  to  entrainment  of  surrounding 
fluid  or  wake  shedding)  signif icantly  in  the  time  it  takes 
the  vortex  ring  to  pass  the  measuring  point.  This  is  justi¬ 
fied  by  flow  visualization  studies  of  the  variation  of  ring 
speed  Uq  and  size  of  the  ring.  The  radius  of  the  ring  r  is 
defined  as  the  point  of  zero  velocity  along  Z  =  0  in  the°steady 
coordinate  system  of  Fig.  14.  Also  note  that  all  the  data  to 
be  presented  are  on  the  "bottom  half"  of  vortex  ring  as  shown 
in  Fig.  14.  A  summary  of  the  data  is  given  in  Table  1. 

Two  typical  data  curves  for  velocities  in  the  vicinity  of 
the  core  of  Ring  No.  1  are  shown  in  Fig. 15.  An  increase  in 
scatter  over  the  data  in  Fig.  13  is  evident  due  to  time  delay, 
jitter,  slight  variations  in  the  ring  position,  and  possibly 
velocity  fluctuations  in  the  core.  Prior  to  data  analysis 
the  Polaroid  oscilloscope  pictures  are  photographically 
enlarged  2.5  times  and  a  smooth  curve  drawn  through  the  numer¬ 
ous  dots.  The  results  of  this  enlargement  are  shown  in  Fig.  16 
and  17.  These  curves  represent  a  portion  of  the  data  taken  on 
Ring  No.  1.  The  remaining  data  on  Ring  No.  1  were  taken  in  the 
vicinity  of  the  core  with  an  expanded  velocity  scale. 

The  U  component  of  velocity  versus  r  along  Z/r  =  0  for 
the  two  vortex  rings  is  shown  in  Fig.  18.  The  corresponding 
curve  for  Hill's  spherical  vortex  is  shown  for  comparison  in 
this  and  the  following  graphs.  The  value  of  rfi  the  ring 
radius  and  the  core  radius  are  defined  by  this  curve.  The 
core  radius  a  is  defined  as  one-half  the  distance  between  the 
positive  and  negative  peak  velocities.  This  gives  a  ratio  of 
core  radius  to  ring  radius  of  a/rQ  =  .27  for  Ring  No.  1  and 

a/rQ  =  .075  for  Ring  No.  2,  indicating  that  Ring  No.  2  is  a 

relatively  thin  core  ring  while  Ring  No.  1  has  a  relatively 
fat  core  but  not  nearly  as  fat  as  Hill's  vortex  where  a/r  = 
1.414.  No  data  were  taken  in  the  core  of  Ring  No.  2  due  to  the 
lack  of  smoke  particles  in  the  vicinity  of  the  core.  This 
defect  of  smoke  particles  is  shown  in  Fig.  12  which  is  a  ring 
very  similar  to  Ring  No.  2. 

A  circulation  T  is  obtained  by  integrating  the  veiocity 
along  the  curve  C  of  Fig.  14.  The  curve  is  a  rectangle  with 
one  side  a  line  along  r  =  constant  and  the  other  sides  tend¬ 
ing  to  infinity.  Thus 

y  -  J ~  (uc?)- Uo)d2. 

O  N 

r=  conn. 

In  order  to  correct  for  the  finite  extent  of  the  data  it  is 
assumed  that  (U(z)-U  )  -1/Z3  after  the  last  data  point,  which 
gives 


17 


mmmmm 


n 


r  = 


j  (uU)-Uo)d(Z  -ita(e,)Zz  ■ 
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where  Z_  is  the  upstream  cutoff  and  Z.  is  the  downstream  cutoff 


of 


the 
2  is 


The  circulation  was  determined  by  finding  the  area  under 
the  curves  of  Fig.  16  using  a  planimeter  and  adding  a  small  cor¬ 
rection  (<i5%)due  to  the  finite  extent  of  the  data.  The  total 
circulation  I*  is  the  integral  along  r  =  0.  The  result 

these  calculations  is  shown  in  fig.  19.  The  meaning  of 
curves  is  similar  to  that  of  Fig.  18  indicating  Ring  No. 
the  thinnest  core  ring  having  the  circulation  more  concentrated 
near  r/r  =1.0.  The  circulation  calculated  ir  this  manner  can 
give  some  indication  of  the  vorticity  distribution  by  differ¬ 
entiating  the  curves  of  Fig.  19.  However,  in  order  to  obtain 
the  vorticity  this  derivative  would  have  to  be  multiplied  by 
an  appropriate  length  corresponding  to  the  distribution  of  vor¬ 
ticity  in  the  Z  direction.  That  is,  the  slope  of  the  curves 
of  Fig.  19  represents  the  vorticity  weighted  by  an  axial  dis¬ 
tribution  length  which  is  unknown. 


Since  complete  two-dimensional  data  were  taken  on  Ring 
No.  1,  further  data  analysis  is  possible.  The  streamlines  and 
vorticity  distribution  throughout  Ring  No.  1  are  calculated 
be  1 ow . 


The  streamlines  are  found  by  plotting  lines  of  constant 
streamf unction  0.  The  streamf unction  C'  is  calculated  by 
integrating  the  velocity  along  appropriate  curves.  For  ar. 
axi symmetric  flow 


Integrating  from  Z  =  +®  along  a  line  r  =  constant  gives 

t! 

(ft  (  r«  copjt,  r.  ocij  =  -  Y'  [  v&)JLl' 


dll-  =  rUdr  -  rVdz 


Assuming  that  V  vl/Z3  after  some  cutoff  Z.  and  non-dimens ional- 
izing 


-  z-\ v i _  f2  y(%)' 

7  ua  tlrJ  n,  L  zu.rc  J3  u.  r°_ 


(2) 


The  s treamf unction  is  thus  determined  by  integrating  the  curves 
of  Fig.  17  with  the  results  shown  in  Fig. 20. 


The  s treamf unction  along  Z  =  0  is  found  by  integrating  the 
U  component  of  velocity  in  Fig.  18. 
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The  results  of  this  integration  are  shown  in  Fig.  21. 

shown  are  the  values  of  - — ^ calculated  from  Eq .  2. 

r 

o  o 

critical  check  shows  that  the  overall  consistency  of 
is  excellent. 


Also 

This 

the  data 


The  actual  streamlines  are  shown  superposed  on  a  photograph 
of  Ring  No.  1  in  Fig.  22. 


The  vorticity,  ,  is  calculated  along  Z  =  0  by  differen¬ 
tiating  the  data 


co  c  »\  o') 

( “vj 


9 

_ 

i.0  5<-SW 

The  first  term  i?  the  slope  of  the  curves  of  Fig.  17  at  Z  =  0 
and  the  second  term  is  obtained  by  differentiating  the  smoothed 
curve  of  Fig.  18.  The  results  of  differentiating  (Fig.  23)  show 
scatter  in  the  vicinity  of  the  core.  This  is  caused  in  part 
by  the  directional  ambiguity  of  the  LDV  system.  The  V  compon¬ 
ent  of  velocity  rapidly  changes  from  a  large  positive  value  to 
a  large  negative  value  giving  data  curves  (Fig.  15b)  in  which 
the  determination  of  the  slope  at  V  =  0  is  difficult. 


The  vorticity  distribution  (Fig.  24)  shows  the  expected 
concentration  of  vorticity  in  the  core.  The  vorticity  at  any 
point  in  the  flow  field  is  known,  since  w/r  =  F(^),  in  a 
steady  axisymmetric  flow.  The  function  F(^)  is  shown  in  Fig. 
25.  Again  this  shows  scatter  in  the  values  of  u/r  for  r/rQ  <  l 

and  r/i'o  >  1;  however,  the  shape  of  the  resulting  curve  is 

the  same  in  both  cases. 


The  problem  of  stability  of  vortex  rings  was  treated  in 
a  separate  study  and  reported  in  Ref.  6.  The  conclusion  of 
Ref.  6  is  that  vortex  rings  are  unstable  to  azimuthal  pertur¬ 
bations  The  measurements  presented  here  were  all  taken 
prior  to  instability. 
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CONSTRUCTION  OF  A  ROTOR  FOR  USE  WITH  LDV  SYSTEM 


A  two  bladed  rotor  with  hollow  rotor  blades  was  designed 
and  constructed.  The  hollow  blade  permits  smoke  to  be  deposited 
at  the  tips  or  any  other  place  along  the  blade. 

A  cross  section  of  the  blade  in  Fig.  26  shows  the  epoxy¬ 
fiberglass  structure  with  a  stainless  steel  spar,  beryllium- 
copper  inner  chamber  and  a  tungsten  rod  to  maintain  balance 
about  the  quarter  chord.  Initially  an  aluminum  model  was 
machined  from  which  a  mold  (Fig.  27a)  was  cast  of  Stycast  2850, 
a  high  strength  casting  resin. 

The  beryllium-copper  was  spot-welded  to  the  stainless  steel 
spar  and  the  tungsten  rod  silver  soldered  on  the  front  of  it. 

An  inner  mandrel  was  inserted  inside  the  beryllium-copper  and 
the  outside  wrapped  with  fiberglass  wet  with  Eccomold  L-28 
laminating  resin.  This  structure  was  compressed  in  the  mold  to 
form  the  NACA-0012  airfoil  section.  After  the  epoxy  cured,  the 
blade  was  removed  from  the  mold  and  the  inner  mandrel  removed. 
The  resulting  hollow  blades  are  shown  in  Fig.  27b  and  27c. 

The  drive  stand  for  the  rotor  (Fig.  28)  has  a  hollow  shaft 
to  allow  smoke  to  pass  to  the  rotor.  The  shaft  is  driven  by  a 
1/2  h.p.,  10,000  rpm  motor. 

At  the  present  time  the  rotor  has  been  operated  at  speeds 
up  to  6,000  rpm.  Velocity  measurements  in  the  wake  of  this 
rotor  will  be  presented  in  the  next  contracting  period. 
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APPENDIX 


MOTION  OF  SPHERICAL  PARTICLES  IN  A  FLOW  FIELD 

AN  LDV  measures  the  velocity  of  particles  in  a  fluid  flow. 
In  the  present  experiment  oil  particles  were  generated  by  a 
Farval  vortex  mist  lubricator,  a  device  normally  used  for  lubri¬ 
cating  high  speed  bearings.  The  oil  particles  are  generated 
by  mechanically  breaking  up  oil  drops  at  a  sonic  orifice.  This 
results,  after  filtering,  in  cold  oil  particles  in  the  micron 
diameter  size  range.  This  method  of  obtaining  particles  was 
found  to  be  extremely  convenient. 

In  order  to  estimate  the  difference  between  fluid  and  par¬ 
ticle  velocities  the  equation  of  motion  of  the  particle  is  non¬ 
dimens  ionalized  and  the  appropriate  dimensionless  parameters 
derived. 


The  equation  of  motion  of  a  small  spherical  particle  sub¬ 
ject  to  Stokes  drag  is  (Ref.  7) 

^  d3  fp  =  -3TT-V)* 


(  A  flp  _  cL  USr  ^ 
^  dt  dt  J 


where  U 
o 
v 
d 


=  velocity 
=  density 
=  viscosity 

=  diameter  of  particle 


-  IT  JL2  fr 
z  c,  r* 


Subscripts  p  and  f  refer  to  particle  and 
fluid  respectively. 


Non-dimens ionali zing  with  respect  to  a  reference  fluid  velocity 
U^p  and  length  L^pp  and  using 


Ap  =  -  p 


d^ 

f  dt 


St  k  ^p  B 


( A-  1 ) 
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where  all  quantities  are  nondimensional  and 


^npp  P f 

Stk  =  Stokes  number  =  - -  (1  +  — ) t 

lref  V 


t"  = 


_  cCZ 


/Si) 


•  (— ) 


In  the  present  study  oil  particles  are  used  in  air,  so 


implying  B  =  0  and  Stk  = 


The  role  of  the  time  scale  t  can  be  seen  by  considering 
the  simple  problem  of  a  particle  moving  in  a  straight  line 


with  a  velocity  at  t  =  0  of  Uq  in  a  stationary  fluid.  At  any 


time  t  later  the  velocity  will  be 


«L-. 

U0 


showing  that  t  is  the  characteristic  time  for  a  particle  to 
reach  the  local  fluid  velocity. 


For  one  micron  diameter  oil  particle  in  air  t  =  4  x  10  sec 

and  the  Stokes  number,  which  is  a  measure  of  the  difference 
between  fluid  and  particle  velocities,  is 


Stk  =  \  x  10  for  U  =  1C  ft/sec  and  L  =  .01  ft. 
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APPENDIX 

MOTION  OF  SPHERICAL  PARTICLES  IN  A  FLOW  FIELD 

AN  LDV  measures  the  velocity  of  particles  in  a  fluid  flow. 
In  the  present  experiment  oil  particle >  were  generated  by  a 
Farval  vortex  mist  lubricator,  a  device  normally  used  for  lubri¬ 
cating  high  speed  bearings.  The  oil  particles  arc  generated 
by  mechanically  breaking  up  oil  drops  at  a  sonic  orifice.  This 
results,  after  filtering,  in  cold  oil  particles  in  the  micron 
diameter  size  range.  This  method  of  obtaining  particles  was 
found  to  be  extremely  convenient. 

In  order  to  estimate  the  difference  between  fluid  and  par¬ 
ticle  velocities  the  equation  of  motion  of  the  particle  is  non- 
dimensio.ialized  and  the  appropriate  dimensionless  parameters 
derived. 


The  equation  of  motion  of  a  small  spherical  particle  sub¬ 
ject  to  Stokes  drag  is  (Ref.  7) 

"b  jS  ^  3?  * 


-  TT  JL*  ft 
z  6  r* 


(  ^  _  d  ^ 

^  dt  dt  J 


where  0 

P 


v 

d 


velocity 
dens i ty 
viscosity 

diameter  of  particle 


Subscripts  p  and  f  refer  to  particle  and 
fluid  respectively. 


Non-dimensionaiizir.g  with  respect  to  a  reference  fluid  velocity 
and  length  L^p  and  using 


Lp  =  - 


Slk  ”  B  +  (Up  O  (a- i ) 
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where  all  quantities  are  nondimensional  and 


Urfr  pf 

Stk  =  Stokes  number  =  - -  (1  +  — )  t 

LREF  Pd 


t-  -  —  I  Pf 


-  /8  A  (Tr) 

In  the  present  study  oil  particles  are  used  in  air,  so  — 


implying  D  =  O  and  Stk  = 


The  role  of  the  time  scale  :  can  be  seen  by  considering 
the  simple  problem  of  a  particle  moving  in  a  straight  line 
with  a  velocity  at  t  D  0  of  Uq  in  a  stationary  fluid.  At  any 

time  t  later  the  velocity  will  be 

Up 

tfo 

showing  that  t  is  the  characteristic  time  for  a  particle  to 
reach  the  local  fluid  velocity. 

For  one  micron  diameter  oil  particle  in  air  -  -  4  x  lO-6  sec 
and  the  Stokes  number,  which  is  a  measure  of  the  difference 
between  fluid  and  particle  velocities,  is 

Stk  =  4  x  1C  3  for  UR£F  =  10  ft/sec  and  L^p  =  .01  ft. 
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Table  1 


DATA  SUMMARY 


6.07  3.10  37900 


The  core  diameter,  2a,  is  defined  as  the  distance  between  the 


/  curve  (Fig.  18) 
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Doppler  shift  vs.  velocity  for  lie-Ne  laser  at  various 
scattering  angles. 
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Fig.  4a.  Doppler 
signal  from  .0002 
dia.  wire  passing 
through  center  of 
probe  volume. 

Vert. scale  lOOMV/cm 
Horiz.  "  .lMS/cm 


Fig.  4b.  Doppler 
signal  from  .0002 
dia.  wire  which  did 
not  pass  through 
center  of  probe  vol. 
Scale  same  as  4a. 


Fig.  4c.  Doppler 
signal  from  spin¬ 
ning  disk. 

Vert. scale  lOOMV/cm 
Horiz.  "  50MS/cm 
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Fig.  5.  Signals  from  Ref.  Beam  &  Dual  Scatter  Systems 
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ISmW  He  Ne  LASER 

laser! - -7 


PRISM  TYPE  \ 
BEAM  SPUTTER 


LOUD  SPEAKER 


PERIODICALLY  PRODUCED 
VORTEX  RINGS 


MEASURING  POINT 
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PHOTO¬ 
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PINHOLE  APERTURE 
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AMPLIFIER 


FUNCTION 

GENERATOR 


SYNC  (Fiq.9o) 
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VARIABLE 
TIME  DELAY 


SAWTOOTH 

SWEEP 

GENERATOR 


(Fig.  9c) 


Tig.  6.  One  Dir-onnlonnl 
LDV  Syntora. 


SWEEP  CPFPTRIIM  SIGNAL 
INPUT  .ANALYZER  .  »NPUI  (Flfl,9b) 


OUTPUT 
(Fig.  9d) 


GATE 


TRIGGER 
_ OUT 


STORAGE 

OSCILLOSCOPE 


Beam  pattern  for  two  dimen¬ 
sional  LDV  looking  toward 
laser.  The  arrows  indicate 
the  direction  of  polariza¬ 
tion.  Beams  1  and  2  measure 
velocity  component  in  the 
x  direction  and  beams  3  and 
4  in  the  y  direction. 


GATE 


Fig.  10.  Demodulation  of  Test  Signal 

f  =  500  kHz 
o 

Df  =  450  khz 
f  =  33  cy/sec 
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DISTANCE  ,  x  (INCHES) 


19.75  FT/SEC 


19. 1 6  FT/SEC 


2.987c 


Linearity  and  Velocity  Chock 
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U  component  of  velocity 

at  r/r  =  .876  for 
o 

Ring  No.  1. 

Time  Scale  20  MS/cm 


Fig.  15  b.  V  component  of  velocity 
at  r/r  =  -■  1.0  8  for 


Ring  n8.  1. 

Time  Scale  20  MS/cm 


component  versus  axial  distance  at  various  radial  stations 
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RING  *2 


Calculated  from  Eq.  2 


Solid  curves  calculated  from  Eq 


HILL'S  SPHERICAL 
VORTEX 
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Fig.  27b.  Finished  blade 


Fig.  27c.  End  view  of 

finished  blade 
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